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‘cate 2 MeO 


An experimental study of the diffusion of helium in the wake of a 
circular cylinder was conducted in the GALCIT hypersonic winc tunnel 
ata Mach number of 5.3. The cylinder was constructed of material 
having random porosity and was mounted with its axis perpendicular to 
the stream. The light gas was injected in small amounts and the thermal 
conductivity method was utilized to detect the concentration of helium in 
the air at points downstream. Problems in the utilization of the thermal 
conductivity method for low sample densities were overcome by suitable 
calibration. 

Flow in the wake of the cylinder was found to display character- 
istically similar behavior at a few diameters downstream, with respect 
to decay and spread of the concentration. Reynolds number similarity 
was established in the laminar case, but turbulent Reynolds number 
similarity may require reference to momentum thickness, which was not 
possible with the present data. 

Profile data was sornewhat marred by a tunnel pressure perturbae- 
tion, but many of the important conclusions were not affected. The pro- 
files appear to follow the theoretical Gaussian distribution in the similar 
region. 

The thermal conductivity method is quite promising as a means 
of tracing the diffusion of one binary gas constituent in another, as 
applied to hypersonic wind tunnel experiment. it will also serve in the 
analysis of transition and turbulence, and of the lateral spreading of 


the turbulent fluid into the rest of the wake region behind the bow shock. 
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I INTRODUCTION 


The problem of the wake behind a blunt body is one of the oldest 
in the classical, incompressible aerodynamic literature. The practical 
aspects have centered around the effects of wake vorticity and circulation 
and the momentum defect determination of drag**, A limited quantity of 
this classical literature has been directly concerned with the problems 
of mass diffusion in the wake. There exists, however, @ great amount 
of theory and experimental data on jet mixing of incompressible flows 
and lower Mach number compressible tlows. The jet mixing theory is 
of course directly analogous to wake theory in many respects, especially 
in the subsonic case’. Also, for the viscous-temperature effects of 
primary interest here, namely those of the diffusion of momentum, 
mass and energy, the jets and wakes are just free boundary cases of 
boundary layer flow for which there is extensive treatment in the 
livergture’. | 

A summary of the mixing theory and list of references to 1954 
has been contributed by pai. This compilation indicates that the know- 
ledge of diffusion processes of incompressible flow is fairly complete, 
and consists for the most part in consideration of momentum-vorticity 
diffusion. The compressible subsonic theory is also well established, 
with the addition of some topics of heat energy difiusion. More recent 
contributions which are of particular interest in the present investigation 
are the works of Tonnes and Schubauer and TEneee: Townsend 
proposes such ideas as that of the intermittency at outer boundaries of 


turbulent flows and discusses the wake behind a two-din:ensional 





* Superscripts denote references at the end of the text. 
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cylinder at some length. Schubauer and Tchen collect and extend the 
ideas of Townsend and others in a well-integrated description of such 
flows, and include an excellent bibliography of the associated literature. 
An additional recent contribution to the specific field of interest here is 
the book by nee’. 

Literature on the supersonic diffusion processes in wakes cannot 
depend so heavily on the jet mixing literature, and very little has been 
done, even recently, in this more specific field. Also, the literature 
described indicates little concern with mass addition and chemical 
reactions, which have recently become highly important in the hypersonic, 
high altitude flight regimes of advanced weapon design and space explora- 
tion. A number of theoretical and experimental papers have been written 
in recent years which have demonstrated clearly the importance of these 
phy sical-chen.ical considerations in the iransfer of heat to bodies in 
high Mach number flow mete A good summary of associated work 
is given by nee. Rubesin - and Wenner. , in particular, show that 
a light gas is most effective with respect to weight of material injected 
in reducing heat transfer rates. One concludes frorn sucha literature 
survey that many hypersonic problems of practical nature will involve 
the injection of materials in light gaseous state into the boundary layer 
for cooling, which in turn may result in chemical reactions and nonuni- 
form mixtures of gases in the wake flow behind the body. Any complete 
study of the hypérsonic wake must consicer these phenornena. 

The wake flow behind a hypersonic body involves the usual cone 
siderations regarding afterbodies, structure, and control surfaces, but 
it hae other implications as well which may become quite important. 


There is, for instance, the probable radar reflectivity and optical 





emission properties of the material in the wake, which are of interest 

in the study of meteor wakes and in the tracking of reentry vehicles, in 
addition to other obvious military considerations. (Dr. Feldman and 

his associates at AVCO Research Laboratories are currently investiga- 
ting these problems.) On the other hand, the study of the hypersonic wake 
is of current interest as a part of the present exhaustive investigation of 
the hypersonic flow field in general. 

The above considerations have led to the establishment at GALCIT 
of a program for systematic study of the hypersonic wake behind blunt 
bodies. Wake pressure measurements at a nominal Mach number of 
5.3 have been undertaken by J. M. McCarthy, and temperature and hot- 
wire anemomeitry measurements are being conducted by A. Demetriades 
and C. F. Dewey. Ina separate study, M. D. Cofiin is continuing the 
work of McMahon” in the investigation of heat transfer with mass 
injection at the stagnation point of a blunt body. In his work, Coffin 
has developed an apparatus for the analysis of the concentration of one 
gas in another by the method of thermal conductivity. The method was 
known to him through his contact with the work of Rush and For cial 
and Forstall and Shapiro: ’, at the Massachusetts Institute of Technology. 

The present investigation was suggested by Captain Coffin (USAF) 
and Professor Lees, as a means of extending the GALCIT wake program 
by a diffusion study. The use of a tracer gas should help in understanding 
the chemical aspects of wake flow and difiusion in a binary mixture 
of gases. It was also envisioned that a light pas might be usedaga 
tracer to determine turbulent processes in the wake. Although the 
thermal conductivity method is a classical one in gas analysis’ °, it 


has seen very little employment in the field of wind tunnel experiment. 





In view of some of the limitations of other devices in the usual hypersonic 
wind tunnel apparatus for the measurement of the physical state of the 
gas, it has become evident that new methods for obtaining and correlating 
this type of data are required. <A large portion of the discussion below 
will necessarily involve the setting down of simple procedures, 
corrections and considerations that have been found essential in the 
employment of the apparatus. It is hoped that inclusion of all relevant 
procedures will assist future investigators in avoiding the many small 
technical delays that interfere with a smoothly conducted experiment. 
Some improvements to the equipment which could not be incorporated 
in the present study will also be indicated. 

Thus, the present investigation is a study of the effectiveness 
of the thermal conductivity method for determining concentration of a 
tracer gas. At the same time it is an exploratory study of the utility 
of this tracer gas in uncovering the nature of hypersonic wakes behind 


blunt bodies. 








II. EQUIPMENT AND PROCEDURE 
li.1. General Experimental Method 


In accordance with the cited objectives the following experimental 
steps were planned: 

(1) Ejection of light gas from a porous body of simple geometric 
shape at a nominal Mach number of 5.8 in carefully controlled small 
quantities. 

(2) The extraction of samples of the wake at various stations in 
the tunnel by means of a small pitot tube. These samples were to be led 
to a thermal conductivity cell which would compare concentration of the 
diffusing agent in the sample with a known reference gas. This procedure 
requires a calibration of the thermal conductivity of known mixtures of 
the two gases against cell readings. 

(3) The correction of the data for mean tunnel concentration, 
differences in pressure of the sample and reference cell, instrument 
errors, and any other effect that might alter consistency and reproduci- 
bility of the readings from test to test. 

(4) The comparison of wake pressure data taken for the ejection 
model against that of the solid model used for total head surveys, 
(McCarthy's model), to correlate results. 

(5) Experimental check of the two-dimensionality of the tunnel 
flow at the center of the model, and spanwise uniformity of the tracer 
gas ejection. 

(6) The reduction of the data to find what conditions of similarity 
may be found in the flow, and whether turbulent zones can in fact be 


defined by the method. 





(7) The comperison with available theory and related experiment. 
ll.2. Selection of Diffusion Gas 


The selection of a tracer gas for diffusion into the tunnel air 
that would best meet the requirements of the experizrent was simple. 
If one checks the list of the common unobjectionable, relatively 
inexpensive laboratory gases, helium and hydrogen are the two most 
suited to the application and most sensitive to thermal conductivity measure- 
ment. Helium is 5.97 and hydrogen 7.15 times ag conductive as air, but 
hydrogen is eee aie: (See Table I.) The choice of helium becomes 
obvious, and it is certainly one of the light materials that rnight be 
employed in surface cooling applications. The commercial helium 
utilized was 99.9976 per cent pure, with traces of carbon dioxide, 


argon, hydrogen, nitrogen, and rmethane, by typical analysis. 


I.3. De scription of the Wind Tunnel 


The GALCIT hypersonic wind tunnel is a closed return, contin- 
uously operating tunnel with two test section legs. The legs are designed 
to be used alternately to provide for installation in one while the other 
ig operating. Leg 1 has nozzle blocks set for a nominal Mach number 
of 5.8 in the 5'' x 5.25"! x 29"' test rhombus. The leo 1 operating 
limits are as follows: reservoir pressure, P= -5 to 100 psig; 
reservoir teniperature, 26 = 225 to 3Z5°F; for which Reynolds number 
per inch, Re = 3.88 x 107 to 30.3 x 10° - The Reynolds nurnber per 


unit length is obtained by the formula 





a ma Sen 
Re/L =(Po/zo) M (7-1) CLT (T Ja 1 , (la) 


O 
for a perfect gas, where zero subscripts refer to reservoir guantities, 
iseatropic expansion from PG ig assumed, and the cquation for the 


variation of viscosity with temperature ig the Sutherland formula: 


ie mre (1/7997? ) (2) 


For these tests with M= 5.8 and - eS ae 


Refin = 2260 x RP. ; (1b) 
where BP is measured in psia. The velocity of the test section flow at 
this nominal Mach number and T. is U = 2776 it. /sec. as computed by 
the formula: 


3 ; 
(U/a,)* = = a eee = NIRT, ; (3) 
(2/M") + (¥- 1) 
where a, is the reservoir speed of sour“. The reservoir pressure 


could be varied between approximately 0 - 100 psig, but critical 


Se 
starting and running conditions were avoided by limitations of 10 - 85 
psig. 

A more complete description of the tunnel and compressor plant 
is to be found in References 19 and 20. A schematic drawing of the 


tunnel is included in Figure 1 and of the test section and related 


experimental equipment in Figure Z. 





ll.4. Model Cons:ruction and Installation 


The basic model of the GALCIT hypersonic wake study at present 
is the two-dimensional circular cylinder. The body geometry was 
purposefully kept as elemeniary as possible in order to simplify the 
correlation of data. Once the complete flow field for this model is 
established the experimental n.ethods for the sin.ple model can be 
extended and applied to r1.ore complex shapes. it was found thata 
0.3'' cylinder diameter was the maximum that could be employed for 
consistent starting of the tunnel flow. Thus, a.3" x 5"' cylinder was 
mounted horizontally frem wali to wall in the forward part of the 5 inch 
wide tunnel test section (Figure ¢). Initially the n:odel was installed at 
a point 13.2'' aft of the nozzle thioat. In the later tests, however, the 
m:odel was moved an additional 2.6" aft to obtain betier flow conditions. 


(See the discussion in Section III. ) 





McCarthy's total head and static pressure surveys were being 
made witha brass model. Although brass might be made porous, ii 
was thoughi that a more suitable n.aterial should be found witha surface 
as smooth as possible anc with natural random porosity. These proper- 
ties were found in an alumina reiractory rod of ihe type used to support 
furnace heater elements. The refractory rods initially were of . 3&5 inch 
exterior diameter and 0.15 inch inierior diameter. The exterior of the 
model was carefully turned toa sn.ooth, uniform .3 inch diameter, but 
the interior diameter could be as muchas. 0Q0Q3 inch non-concentric, 
because of warping of the rod. The slight non-symmetry of helium flow 
that would result from this interior eccentricity was considered of nepli- 


gible concern as to its effect on wake measurements of gas concentration, 





in view of other experin_ ental difficulties that appeared to exclude such 
minor considerations at the time of model selection. Unfortunately, the 
rod material was somewhat brittle and the models had to be replaced 
frequently because of breakage. In spite of these deficiencies of the 
alumina rod the material was considered well suited to the present 
application, because it provided low, well distributed helium flow rates 
with low metering pressure during bench tests. Thus,ihe momentum of 
the ejected gas is small for a given mass flow, and should not appreci- 
ably alter the external flow and shock wave structure. It was desirable 
to achieve this result, if possible, in order to correlate diffusion data 
with the pressure data of McCarthy's experiments as directly as 
possible, and therefore obviate the necessity of duplicating a large 
amount cf his work. 

The bench tests of the models mentioned above included other 
considerations which are peculiar to this experiment. First, it was 
noted during preliminary tunnel runs that a lightesurface coating of 
oil would collect on the model surface. The oil problem is a continuing 
difficulty in the GALCIT tunnel for which a satisfactory solution could 
not yet be found. It was therefore necessary to measure in sore way 
the effect of the oil on the flow of helium through the rod surface. After 
several hours of running, ihe lightly oil-coated side of the rod was 
covered and the model was connected to a vacuum pump at one end, with 
the other end sealed. The difference in the amount of vacuum the pump 
would draw against a sealed space, and against the uncovered side of the 
model, was measured with a U-tube manometer. The oil coated side 
showed a negligible change in porosity. The oil problem therefore was 


not considered serious, since the model could be rotated from test to 
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test to collect an even distribution of cil, and the flow level could be 
adjusted to provide the same flow of heliwyr. regardless of slight 
changes in model porosity. 

One 1..Uu5t also rake certain that all the metered hclium is 
ejected throuch the cylinder and aot lost in faulty seals at the tunnel 
walls. A clamping and locking arrangement and soft 'o''-rings were 
installed to prevent such leakace. After each installation the vacuum 
test described above was applied to determine if the seal was con- 
sistently good. The helium was injected from both ends of the model in 
order to provide a more uniform flow distribution across the cylinder 
span. The rnodel was mounted between metal ports because of the 
difficulty of securing an injection model in glass ports when the in- 
stallation must be repeatedly disassembled. Schlieren studies were 
therefore not teasible with the cylinder injecting helium. However, 
comparison of the flow characteristics with and without helium flow 


could be obtained by taxing representative total heac pressure traverses. 


1.5. Control of Heliurn Injection 





The ilow of helium entering the cylinder at the tunnel walls was 
metered witha standard risher-Porter Tri-flat flowrator, (Tube No. 
3-F-3/8-25-5/70, having a 3/5 inch glass ball), in series combination 
witha smaller Fisher-Porter flowmeter (Tube No. OIN-15, 1/8 inch 
steel ball). Zhe large tlowrator provided stable measurement of the flow 
in spite of "back-pressure'' variances of the model, while the small 
flowmeter provided a double check and a rmore sensitive, vernier-type 
monitoring of the flow quantity. The metering pressure was maintained 


by reference toa standard U-tube mercury manometer. (See Figure 2.) 





1] 


Although the temperature correction is smell, the temperature of meter- 
ing was noted for each test. 

The metering of helium through the flowmeters was checked 
against the Fisher-Porter Company predictions and the previous 
labcratory calibrations for helium: flow. This calibration was accome- 
plished by the timed cisplacemeni of a large volume of water at several 
heliurn flow rates. The large flowrator checked closely with the 
manufacturer's prediction, but the smaller flowrator was affected too 
seriously by back-pressure of the water to accomplish a good test. 

The flow of helium was therefore always based on the reading of the 
large flowmeter, with the smaller meter as a check against variances 

in cylinder porosity by the back pressure cifect. The test was conclusive 
to about = 3.5 per cent accuracy in the flow range of interest, whichis 
also approximately the crder of accuracy in setting of the flow for the 


tunnel experiments. 


Il.6. Sampling Probe 


The probe employed to withdraw flow samiples from the wake of 
the injection model was of the type designed for total pressure measure-= 
ments in boundary layers. The tip was constructed by compressing a 
hypoderniic tube to a flat orifice of 0. 004'' depih with 0.039'' width. The 
tip was faired out to a 0.25'' ©. D. stainless steel tube fitted witha 
0.067'' I. D. neoprene tubing at the downstream end for transport to the 
gas analysis cell. The same probe was used for total head traverses in 
the wake. 

The probe was mounted in a mechanism designed to traverse it 


axially from zero to 27 model diameters aft of the model centerline, and 
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Sen aed 


vertically 1. 2"' upward io 2.0" below the axial centerline of the tunnel. 
The probe could be set to an accuracy of about 0.01" axially and 0. 001" 
vertically. (Sce Figure 2.) 

Samples were withdrawn despite the low prescure of tne tunnel 
by the creation of a near vacuum condition in the concuctivity cell. 
(This procedure will be cxplained in more detail in the next section. ) 
In the case of the lower tunnel pressures (near the model and at lower 
~eservoir pressure) a relatively slow time constant for 
sample was noted - on the order of 30 seconds. It was considered 
essential to utilize the sriall probe, however, because of the narrow 


regions of helium distribution to ve encountered near the model. 


U.7. Thermal Conductivity Apparatus anc Procedures 


The thermal conductivity or heat transfer capacity of a gas mixture 
is a distinct physical quantily that may be used in various ways to detect, 
identify and trace the concentration of one binary constituent in another. 

A very comprehensive summary of the principles and some of the more 
~—— : = 1S pis. ae , 
common applications ig presented by Daynes -. An application in which 
ae | ee es ma a . 16,17 
helium was traced is the work of Rush, Forstall, and Shapire : 
where the method was utilized to trace diffusion in coaxial gas jets at 
low speed. The form of the equipment in the present experiment is 
similar to that designed by Rush and Forstall. It was designed by 
M. D. Coffin at the GALCIT hypersonic laboratory. 

The thermal conductivity method of com.paring two gas mixtures 
consists fundarnentally in passing the mixtures in question through 
identical cell chambers with electrical resistance cler.ents forming legs 


of a Wheatstone bridge. (Sce cell diagram, Figure 3, and wiring diagram, 


wv 
coe 
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Figure 4.) The resistance elements in the two chambers have equal 
currents since the resistances are paired in the two legs of each 
chamber. Small variances can be removed by adjustment of the dividing 
resistance at the bridge input. When the gas conductivity in one chamber 
equals that in the other, the wire heat, (I°R), is conducted through the 
walls of the chamber, in equal amounts for the two chambers, anda 
balance is achieved with no voltage across the bridge. The resistance 
of the wires in the two Chambers is not varied other than by the chamber 
temperatures. Also, in the bridge, one holds the current constant and 
essentially measures a change in R, the resistance of sample chamber 
wires. Thus R= R(T) , where T, is the wire temperature. ee is 
dependent on wire heat logs by the energy balance, 


° ; ; “a 
Gy = gt 4 = IR a 


where neglecting radiation and in stagnant gas: 


af = heat loss through the gas = HI 47/8) 


OT = i - TT. , the difference between wire and chamber 
wall temperature 


6 = typical chamber dimension 
qo = heat loss through wire ends = CaT, where C = coefficient 
of wire conductivity times wire area, in appropriate units. 
Then 
ik 

& m 

= is 4 
IRe (=e + CPT, Te) : (4) 


The chamber walls conduct rapidly enough to assume T, = constant. 
(In this case the chamber walls were constructed of brass.) With cone 
stant current Eq. (4) then gives a relation between R, Bs ees , and Ty, : 


The mixture conductivity, however, is a function of concentration of the 
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mixture species and the temperature, Le , of the mixture, by the 
kinetic theory relations (Cf. Reference 2, page 159). The mixture tem- 
perature of the combination of helium and @ir can be expressed as some 
function of oT, or as approximately the mean, i. = (T - T.)/2. One 
can thus replace K in Eq. (4) by the above equation from kinetic 


theory and then Ty by a linear function of R, arriving at 


R= Rt+aR 5 {K S+4K) - (5) 


Thus, each Wheatstone bridge resistance, KR, is a function only of the 
change in concentration of the helium in an air-helium mixture surrounde- 
ing it. In this case, the bridge resistance in the Cell cavity containing 
helium enriched air is reduced by the greater cooling capacity of the 
light gas. The variation of resistance in one chamber will cause a 
voltage unbalance across the bridge which can be reac ona sensitive 
potention.eter. Equation (5) can best be evaluated in terms of this 
potential for purposes of calibration by laboratory measurenient of 
known sample concentrations. (See Section iI. 7.1.) Other schemes of 
operation are feasible, but this rnethod is one of the moat adaptable to 
this experiment and was especially convenient in that a commercially 
manufactured cell could be utilized. (The specifications and manu- 
facturer of the particular cell chosen are given in Figure 3.) 

Although the fundamental procedures outlined above are con- 
ceptually simple, several additional problems arose in procuring 
representative samples from hypersonic flow fields, and in removing 
possible sources of error from the method. These are the considera- 
tions which determine specific design of the apparaius and operating 


procedure, and which will be discussed in succeeding paragraphs. 








They can be classified under the following headings: 

(1) proper calibration of the cell 

(2) adequate handling of the sample 

(3) pressure differences between the two chambers 

(4) temperature ditference effects 

(5) instrument errors and accuracy. 
(An additional item of consideration here might be the handling of 
dissimilar mass flows between the two chambers. However, since 
hypersonic tunnel samples are widely variable in density and large 
errors would be introduced by convective cooling effects, it was decided 
from the outset to stagnate both sample and reference Saas.) The 
above five items will be discussed in turn with regard to how they affected 


the design and/or procedure of the experiment. 
1.7.1. Cell Calibration 


The calibration curve of potential across the bridge 
versus concentration of helium: in a reference gas must be obtained bya 
laboratory method, as shown in conjunction with Eq. (5). The reference 
gas used must be one which remains consistent in heat conductivity after 
calibration for each subsequent filling of the reference chamber. In this 
case, if one were to employ room air in the reference chamber, filling 
the chamber for each period of use, an error may be introduced by 
departure from the calibration curve because of changes in the atmos- 
pheric mixture. The air mixture may change with small percentage 
variations of carbon dioxide and water vapor. The carbon dioxide 
variation will not be significant, but water vapor partial pressure is 


sufficiently variable and water vapor is sufficiently different in thermal 
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conductivity that noticeable error may be introduced if dry air is not 
used. (See Table I.} Since the air in the hypersonic tunnel is well 

dried and otherwise sufficiently uniform in content frorn one test period 
to the next, a sample of heliumefree tunnel air was taken as the reference 
gas before each test. 

Having selected the reference gas as tunnel air, the cell 
calibration should be made with air-helium mixtures against dry air. 

The collection of large quantities of tunnel air is inconvenient, however, 
and the present calibration procedure involved the com:parison of known 
mixtures of helium in dry nitrogen with pure nitrogen. Since the con- 
ductivity of nitrogen is very close to that of air, this procedure involves 
an error of snaller magnitude than that of preparation of the known mixtures 
and plotting and reading of the graph. The magnitude of the error involved 
can be estimated by noting that nitrogen is 0.996 ag conductive as air, 
(Table J). Thus, an error on the order of one part in 250, or .4 per 

cent, is the result of comparing any sample against a reference of 
nitrogen instead of dry air. The difference in concentration obtained 
between comparison ot a sample containing large quantities of nitrogen 
against a nitrogen reference and comparison of a similar air<-helium 
sarmple against air is therefore guite negligible. 

The samples were prepared by allowing a bottle of high pressure 
nitrogen to discharge into a bottle of helium, with a sensitive gaze 
measuring the pressure of the helium bottle accurately before and after 
the addition of the nitrogen. The partial pressures of the two gases in 
the volume then determined the concentration of the mixture. The 
preparation was therefore independent of possible nitrogen leakage 


during the transfer. The method of computing the concentrations is 





based on the equation of state for components of a perfect gas mixture: 


Prie ~  PHe Re Tin a (Mz Q/ Vin Ro/ M5.) an : 
where 

Pre = partial pressure of He, measured before mixing 

Ro = universal gas constant 

Mire = molecular weight of He 

Ve = mixture volume, constant 

Tin = mixture temperature, constant 

m = mass of helium 

He 
Po = Py + Py. = pressure measured after mixing, by 
z Dalton's law. 
Solving for Mie ? 
Pre ae Mee 


tr = ee 2 me Re 
He Ro 4 


Similarly, 
"Ne Vin me 
Ma) = Sa is the mass of nitrogen. 
2 Oo 7m 

Thus 

K e He = Pue mae (6) 

~ “Hoe me Dat eieetepe= Mame “aa 
ne Met ™N, Pu, VN, * PHe “He 


Equation (6) is utilized to obtain the newly mixed concentration. 
A series of samples was prepared, having different concentrations of 


He in N As each sample was prepared, its concentration was read in 


a 
terms of voltage deflection across the thermal conductivity bridge. The 
plot of these calibration points is a smooth curve with little scatter. 


(Figure 5). 
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Temperature and pressure must be equal in the two cell cavities 
for a proper calibration. (Pressure and temperature are discussed more 
completely in the sections to follow.) Also the calibration curve is appli- 
cable only for the voltage and current utilized across the bridge when the 
sample deflections were read. The voltage and current choice depend on 
two conflicting considerations. Calibrations which are more sensitive 
and linear can be obtained for higher voltage-current combinations, as 
shown by the curves of Rushand For gtan?®. On tne other hand, itis 
desirable to keep cell resistance-wires warm by leaving the electric 
current on while drawing the next sample, in order to minimize the time 
required for the cell temperature to reach equilibrium in the measurements. 
However, during the sample collection, the gas density in the chamber 
becomes very low, and hence the voliage-current combination must be 
low so that the wires will not produce heat faster than it can be dissipated. 
(See Sections II.6. and IJ. 7.2.) Even though the current was turned off 
before gathering a sample, the lowest sample densities are within a few 
millimeters of the vacuum pressure, when the current must be on to read 
the bridge deflection. It is therefore necessary to calibrate the cell with 
the lowest satisfactory power to prevent fusing of the hot wires. To 
accomplish the sensitivity adjustment in this case, a variable resistance 
across the potentiometer terminals is utilized. (See Figure 4.) For 
calibration over the low concentration range of interest in the experiment, 
a current of 80 milli-amps with a 3 volt dry-cell circuit was employed, 
and the potentiometer resistance was set so that an essentially linear 
curve was obtained in the 0 - 2 per cent range. (See Figure 5 for equations 


of portions of the curve. ) 





Y.7.2. Sampling Equipment and Procedure 





The method of obtaining a representative gas sample from the 
hypersonic tunnel is somewhat more critical than in other applications 
because of the low pressures encountered (on the order of oae millimeter 
of mercury). An intermittent sampling system is required in order to 
stagnate both sample and reference for proper comparison. 

As previously mentioned, the reference air was drawn from the 
tunnel in the present experiment prior to commencement of the helium 
flow. The bridge was balanced with the dry tunnel air in both sample 
and reference cavities. Samples were withdrawn, after the helium flow 
was commenced, by creating a near vacuum at the outlet of the cell 
(Valve 3, Figure 2). In order to compress the sample as muchas 
possible to compare it properly with the reference air taken at atmos- 
pheric pressure, a hand-operated mercury pump was included between 
the tunnel and conductivity cell. The pumpin this case had a comprese- 
ion ratio of about three to one. The pump incluced aUe-tube pressure 
scale, accurate to about 1 mm of mercury. Referring to Figure 2, 
valves 1, 2, and 3 are sequentially opened and closed with a continuously 
operating vacuum pump so that the sampling steps are accomplished as 
follows: 

(1) evacuation of the cell, mercury pump, and connecting lines 

(2) washing of the new sample through the same volumes 

(3) re-evacuation of the volumes, which establishes a vacuum 

in the cell 

(4) drawing of sample with valve 3 closed, and the mercury 


flask lowered 
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(5) Compressing of the sample with valve 1 and 3 closed. The 
time required for each of the above steps is a function of the sample 
density. A low density sample enters the cell slowly because of low 
pressure differential between cell vacuum and probe total head. Itis 
important to determine the proper time interval for each step. The time 
of taking one sample point is generally 3-4 minutes at beat by this 
procedure, but if itis hurried the sample may not be representative. 

The time for each step is found by taking an extended time at first, and 
then reducing the interval until a discrepancy is noticed in comparison 
with the same sample point taken over a longer time interval. This time 
interval is most critical for the lowest sample pressures. 

The sampling apparatus was thoroughly leak-checked prtor to 
each test run by applying low vacuum to the mercury pump, tubing and 
cell cavities, closing the valves, and noting whether or not the vacuum 
is maintained by the mercury column. Evena slight leakage of room 
air into the system would cause gross contamination of the sample. Also, 
when @ known sample is put in the sample cavity and left overnight, the 
same reading could be read again before the next test period, asa 
simple day to day check. 

Another pertinent remark should be made in regard to the method 
of collecting the reference sample. Since the pressure calibration 
depends on the difference between the cell cavities, one should be careful 
to see that the reference sample is collected at exactly atmospheric 
pressure. This objective may be accomplished by proper sequencing 
of the valves at the time of collection, (i.e., shutting of the entrance 
valve an instant before the exit valve). Once the reference sample is 


collected at a known barometric pressure, it need not be changed again 





during the test period. 
ll.7.3. Pressure Difference Calibration 


The application of the thermal conductivity method of gas analysis 
to hypersonic flow requires some suitable means of correction for 
differences in pressure between the sample and reference cavities of 
the conductivity cell. For slight pressure differences this error is 
neplizible, but when the reference gas is at atmospheric pressure and 
the sarnple at a few millimeters of mercury absolute, the error in 
neglecting the pressure difference may amount to 50 per cent or more 
of the concentration reading taken. Simple kinetic theory predicts that 
the heat conductivity coefficient of a gas is independent of pressure if 
the mean free path of the gas molecules is much less than the chamber 
dimensions, and proportional to pressure in the opposite limiting case. 
The transition from one limiting case to the other is shown experimentally 
by pigibers”. His craph of heat loss versus iinudsen number 
(proportional to pressure) clearly indicates a smocth transition frorm the 
independent region to the proportional region at Knudsen numbers 
corresponding to the low pressures of the hypersonic wind tunnel, if 
the chamber dimensions are such as to permit sucha transition. It 
was evident from the experiment by Bomelburg and from a few preliminary 
tests that the pressure correction to conductivity would be required in 
the present investigation. 

It is not simple, mechanically, to pump a representative sample 
from: the wind tunnel at near vacuum to a pressure of one atmosphere in 
the sample chamber. Most methods of accomplishing this compression 


involve the use of a pumping fluid medium which may contaminate the 
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sample. While it is certainly desirable to investigate such means 
extensively, a simple expedient to circumvent the need for sucha pump 
is to calibrate the error in bridge potential against the corresponding 
pressure difference. This calibration is accomplished by reading the 
bridge deflection of a known concentration samiple at reduced pressures. 
The sample potential is initially checked at one atmosphere and then 2 
portion removed in steps by a vacuum pump uniil very low pressures are 
reached. The new potential across the conductivity bridge is reaa for 
each step. A suitable means of plotting the resulting calibration curve 
for correction of tests ig evident by physical reasoning, supported by 
experiments like that of Peni One is led to the conclusion that 
the conductivity of the gas in each cell cavity, and therefore any unbalance 
of the bridge, is a function of the absolute pressure in the cavity, when 
the pressure is low enough that wall effects begin to become important. 
Thus, if 

p. = pressure in sample cavity, cm. Hg. absolute 


S 


Pp = barometric pressure of the open and of the U-tube column 
S measuring pressure in the sample cavity, cm. Hg., 
a slowly varying function of time 


Py, barometric pressure, p, , of the room at the time when 
r the reference gas is s collected, cm. Hg. 
4p = height of mercury column, cm. 


Then 


P, = Foo 4£P 


is the independent variable in a function which determines the emf error 
due to low pressure. Now p, = p, for any given calibratiou or test 

S © 
where the barometer does not change more than a few mm of Hg., as 


is invariably the case. The emf correction is not detectable on the 
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instruments when the sample and reference cells are so nearly equal in 

pressure. One can thus say that P, =P, so far as conductivity in the 
e S 

cell is concerned. To provide for barometer fluctuations, the absolute 


pressure of the sample cavity should be piotted in atmospheres as 


(p/P, ) = (p, -Ap)/p, = 1-(op/p,) , 
8 S S 8 


for the abscissa of the calibration plot. The graph of Figure 6 is the 
result of the calibration. The points fall closely along the same curve 

for the three low concentrations utilized. Tunnel air reference points 
cannot be distinguished from nitrogen versus helium-nitrogen points in 
the small scatter at higher pressure differences. The scatter is the 
result of instrurment inaccuracies, sampling variances, and the pressure 
measurement inaccuracy of the mercury pump scale. The pressure 
measurement was accurate to about 1 mm of mercury. Redesign and 
construction for better accuracy was not considered essential in this 
initial wake survey, in view of other experimental inaccuracies of the 
samme order of magnitude. Also, only a few points near the model 

involve pressure calibration at the very low pressures. A maximum of 
ten per cent inaccuracy for the pressure calibration at very low pressures 
was accepted in that if would normally involve much less than ten per cent 
error in the total readings of bridge potential. The corresponding per- 
centage error for any concentration level may be obtained from Figure 6. 
The scatter in the lower concentration region of the curve is entirely 

the result of slight variances in the sample and instruments used at the 
time of the calibration. The emf correction here is small, and the best 
average of the points was taken in the curve fairing. (Later, the 


potentiometer was replaced and the micro-ammeter improved in order 





tJ 
i> 


that leas scatter would oe obtained in taking the wake concentration 
surveys. See the instrument discussion in the next section. ) 

The pressure calibration certainly would show a variation with 
concentration if the concentrations were of a higher magnitude taan the 
3 + 5 per cent levels measured in this investigation. A series of pressure 
difference curves with various concentration levels as a parameter 
would be required for high concentrations. 

Ie is well to note that a pressure calibration for one geometry of 
the apparatus might not be the same as for another where the dimensions 
of the tubing and cavity were changed. This would be the resuit of 
mean free path and cell hot wire end loss cffects. (These effects are 


discussed in the literature of kinetic theory and hot wire anernometry, 


7 


f, Reference 22.) 


Finally, the following list of steps is given for a carefui and 
precise pressure correction of the type that would be required for emf 
accuracies of 97 - 95 per cent with the lowest sample pressures to be 
found near tne r:odel at M= 5.6: 

(1) measurement of sample pressures in the range of 1 - 20 mm 
of mercury absolute to an accuracy of 0.i mm, and higher 
pressures to 1 mm 

(2) accurate calibration curves for known sample concentrations 

(3) collection of reference gas at a known (atmospheric) pressure 

(4) monitoring of barometric pressure throughout the test 

(5) computation of pressure differences between cell cavities 
and normalization by the barometric pressure. {Since the 


variation in barornetric pressure may be as muchas 3 - 4mm 


of mercury for a test run of several hours, each critical 
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reading can be computed using a curve of barometric pressure 
for the times of the test. ) 
The alternative to this procedure is the pumping of each sample 
to atmospheric pressure in the sample cavity with no contamination 
during the process. In the present test, the above five steps were applied 
as Closely as possible, within the 1 mm accuracy of the mercury pump 


pressure scale, and with special care during the critical readings. 


Id. 7. 4. Temperature Difference Effects 


Normal variances in-room temperature will not affect the 
accuracy oi the thermal conductivity bridge as long as the temperature 
is equalized between the sample and reference cell cavities. (This was 
shown by Eq. (5) above, and its derivation.) This statement can be made 
as long as the temperature cocfficient of conduction for the gases to be 
compared is approximately equal. This factor determines the variation 
of heat conductivity with temperature, and is nearly the same for the 
gases considered here (Table I}. The protection from radiating devices 
and stray air currents in the room by a wooden container and the use of 
high heat conductive brass for the cell material, help to insure the 
equalization of temperature. The brass also provides some equaliza- 
tion between cell cavities for samples brought in at a temperature 
unequal to the room temperature of the reference cell. The sample 
traverses a distance of 3 or ¢ feet in the. room and is compressed in the 
mercury pump; also, there is a finite time between sarnple collection 
and reading of the bridge. The temperature of the average sample 
therefore cannot be greatly different from the room temperature of the 


reference cell. It ig well, however, to estimate the effect of a smail 
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temperature difference in the bridge deflection reading. This estimate 
has been accomplished by Davee » who states that the simple pre- 
cautions given above are quite adequate to prevent serious effect of 
temperature differences between cell cavities. One may, however, 
install the cell in a thermostatically controlled housing if very accurate 
readings of low conductivity differentials are desired. 

Daynes, in his book, has givena very extensive coverage of the 
more obscure Considerations in the method of thermal conductivity gas 
analysis. He includes such factors as vibration and other secondary 
effects. His conclusions regarding the detection of highly conductive 
gases in air, however, are that very small concentrations can be 
measured with accuracy if only the precautionary measures of the 


present experiment are employed. 
Wi. 7.5. Control of Instrument Errors 


Referring to Figure 2 the instruments employed for measurement 
of the concentration potential difference were 
(1) the milli-ammeter 
(2) the thermal conductivity cell 
(3) the precision potentiometer with its associated standard 
cell and galvanometer (The galvanometer is utilized to 
obtain a null current through the potentiometer circuit. ). 
A simple method for control of most of the possible errors in these in- 
struments was available through the employment of the samples prepared 
for the calibration procedure. A sample was put in the sample cavity 
before each testi run and the bridge emf was checked. The ratio of the 


calibrated value to the present emf was then a corrective factor for the 
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readings of the bridge during the test. This procedure also served to 
correct for variations in room temperature, as a check on the bias 
setting of the bridge sensitivity, and as a check ior damage or con- 
tamination of the conductivity cell. Along with a standard cell comparison, 
it also served to check the performance of the potentionicter. Carefully 
mixed samples of known concentration are therefore indispensable in the 
use of the thermal conductivity method. Frequently, as meitioned in the 
sample leakage control discussion above, the known sample was left in 
the cell for extended periods of time and the same reading was obtained 
upon recheck. The reproducibility of such readings on known samples 
was within about one per cent at all times. 

The bridge was balanced by a null indication on the galvanometer 
and potentiometer before each test with the dividing resistance at the 
bridge input, (B+, Figure 4). This balance was accomplished with 
nitrogen in reference and sample cavities, and then with air in bothasa 
recheck while gathering the tunnel reference air before injection of 
helium. The stability of the cell zero potential was within 0.05 m.v. 


at all times. 


U.8. Total Pressure Apparatus and Procedures 


The magnitudes of injected helium concentration to be measured 
in a flow depend inversely on the local density and velocity of the flow. 
The determination of these quantities at any station cf interest is 
therefore mandatory to a complete understanding of the diffusicn 
processes. [tis not an easy matter to obtain such quantities ina 
compressible figaa The total and static pressures are relatively 


easy to obtain, but the measuren.ent of another state variable is 





reauired to determine density and velocity. The accitional measurements 
might be the speed of sound, the temperature, or the direct measurement 
of velocity and density, if these measurements were feasible and 
dependable. At present, methods for close determination of these 
flow variables are not fully reliable. One may assume isoenergetic 
flow, however, without serious error in first order calculations, and 
compute an approximate local temperature. The pressures are then 
measured with suitably designed probes and associated manometers 
or transducers, after which the Mach aumber can be obtained through 
the Rayleigh pitot relmion-* 
2¥ 2 YV-l .1/%-1) 
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static pressure 
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total pressure behind probe norrmal shock wave 
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Mach number ahead of probe shock wave 


= 
( 


¥ = ratio of specific heats. 


The Mach number and pressures may then be employed to determine 
other fluid quantities. 

As previously mentioned, a pressure study by J. McCarthy in 
the wake of a 0.3 inch brass cylinder was also in progress at the time 
of the present experiment. It was desirable in terms of wind tunnel 
running time to utilize the pressure data fron: the brass model in the 
diffusion studies made with the porous model. ‘his procedure required 


& comparison of the two models as to flow characteristics, without 
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injection and with helium injection at the level to be metered during 
the wake diffusion tests. The total pressure profiles were considered 
the best means of comparison. Total pressure profiles behind the 
injection model were therefore made with the equipment employed by 
McCerthy. This equipment will be described more completely in his 
report, but a brief description will be given here. 

The probe utilized in the total head traverses behind the cylinders 
was the same instrument employed te gather diffusion samples, and was 
described in Section I].6. The total pressure in the probe was first read 
with a mercury micro-manor.ster at one or two points in the flow, after 
the tunnel had been stabilized at the test reservoir conditions. The total 
pressure was then diverted toa Statham | psi differential pressure trans- 
ducer with a reference pressure provided by a silicon oil column. The 
transducer converted the pressure to electrical potential which was 
traced on a Mosely xy-~ plotter. The transducer-plotter system was 
calibrated by measuring the pressure with the mercury micromancmeter 
ata few points. The output from a voltage divider coupled to the probe 
traversing mechanism was fed into the x- scale of the plotter. Sufficient 
additional micro-manometer points were read to establish that the 
recorder scale had not shifted during the tunnel traverse. All pressures 
were referred to 2 very low vacuum and therefore were plotted as absolute 
pressures. The micromanometer pressures were non-dimensionalized 
with the tunnel reservoir pressure converted to absolute units before 
making the traverse. Thus, the pressure traverses are plotted ina 


convenient form for comparison of flow characteristics. 





30 


lWicCarthy nad also taken static preasure traverses of the wake 
behind the brass cylinder. A comparison of static pressures with the 
porous model, however, was considered unnecessary in that total head 
profiles would give sufficient insight into any flow ditferences, wiich 
should noi be large. 

Results of the wake studies Gepend greatly on whether two 
dimensional ilow actually existed in the tunnel across the cylinder span. 
A few representative total head traverses were thereiore made to check 
the two dimensionality of the flow approximately one inch from the tunnel 


centerline. 
Il..9. Correction and Correlation of Data 


The reduction of emf data from the conductivity bridge requires 
a certain sequence of steps according to the discussion concerning tae 
several items of equipment in the preceding sections. SGefore outlining 
the reduction technique, however, it will be necessary to show how the 
build-up of heliur: in the tunnel circuit may affect the data. 

The quantity of helium in the tunnel air can be predictea bya 
brief analysis. Heliwm concentration coes not build up in value indef- 
initely, but approaches an asymptotic ratio of helium mass to tunnel air 
mass which is dependent on the amount of make-up air added to the tunnel 
circuit and the corresponding leaxage rate frorm the circuit. The 
governing equations are analogous to those for water vapor in the tunnel 
air, where the amount of water vapor is controlled by the bypassing of 
circuit air for drying to prevent condengation in the test section. The 


equation for the tunnel heliurm concentration can be derived as follows: 
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ict m,,, = Constant mass/time of He entering the tunnel 


mn = constant mass/time of make-up air entering the tunnel 
M = mass of mixture in the tunnel, 

K = concentration by mass of He in the mixture 

rm = mass of mixture leaving the tunnel with concentration , 


which must equal the total mass entering for steady tunnel 
operation, or 


mm = et ee = a consfant ° 


Therefore, the rate of change of helium concentration is given Ly 


i Mi.” (Mo, + m.) K 
= pete os eA ee aa ot e 


% 
<t 


The solution of this equation with K = 0 att = 0 is ag follows: 





m (riu.. + m_) 
Ks = i <exp(- t ) (3) 
Mise om I M 


Equation (c) demonstrates that the tunnel concentration tends to an 
asymptotic constant value as t becomes sufficiently large. This result 
was confirmed experimentally in preliminary tests. Within the accuracy 
of the thermal conductivity measurement, the time required to reacha 
constant value was on the order of 10 - 15 minutes for the ncrmal quantity 
of make-up air, ma, , utilized to replenish the tunnel circuit. The tunnel 
concentration with normal make-up air and low reservoir pressure is of 
the same order of magnitude as the lower concentration readings. Itis 
therefore desirable that as much as possible of the background concene 
tration be removed. Make-up air was increased to accomplish this 
removal by excess leakage until levels of helium conceniration about 

.06 per cent or less were obtained at the several reservoir pressures. 


Equation (8) shows that the asymptotic constant value of tunnel 





concentration reached is independent of tunnel -S , and is dependent only 
on the rate of helium and air added to the tunnel. (The rapidity in reach- 
ing the asymptotic value, however, is dependent on Po through the tunnel 
cycle mass ofair, M.) This result was also shown to be correct when 
the valves admitting make-up air were left fixed and the tunnel ae was 
changed for some of the tests. When the mass rate of the helium and the 
corresponding tunnel concentration are known, the make-up air rate can 
be computed by Eq. {8). It the make-up air raie is too high the tunnel 
air will not be sufficiently dried, and errors in measurement of helium 
concentration may be introduced through lowering of the air conductivity 
by water vapor. The tunnel concentration may vary during a test if the 
make-up air is changed. Jt was therefore checked at intervals and 
corrections were made as necessary in the data reduction. 

Perhaps the best method of discussing the reduction of concentra- 
tion test data is to follow a typical data sheet (Table II) and reduction sheet 
(Table IiI), The values of concentration in millivolts read for a series of 
distances, y, from the vertical feeds are entered in the reduction 
sheet, along with the corresponding pressure differences, ap, between the 
current atmospheric pressure and the sample cavity pressure. The com- 
putation is started in column 6 with the subtraction of the tunnel background 
concentration from the emf reading. The tunnel concentration sample is 
always compressed to atmospheric pressure, which allows subtraction 
before Apcorrection. The computation is continued with column 7, where 
the emf reading of the conductivity bridge is corrected for instrument 
errors by the ratio of the emf of a known sample to its current reading. 
(The instrument correction is not large, for it can be controlled by the 
variable resistance across the potentiometer terminals.) This ratio 


corrects for variations in current measurement, room temperature, 





and slight change in cell characteristics as discussed in previous 
sections. 

Column 9 of Table IiI scales the «a p vacuum versus barometric 
pressure reading of the rmercury colunin to that barornetric pressure. 
the barometric pressure taken al one timie during the experiment is 
sufficiently accurate unless very low tunnel densities are expected. if 
these densities are such that the ratio A P/ Pi. >.99, it is well to take 
into account shifts in barometric pressure of greater than 1 mn. 
Fortunately such low pressures are found only very close to the model 
in restricted regions, where the concentration ig normally high. The 
percentage error is therefore low in the total emf reading for a relatively 
high error in Semf caused by cell cavity pressure differences. 

The ap/P, value is utilized in Figure 6 te obtain the Semf 
correction for column 10 of Table III. Finally the emf versus concentration 
graph of Figure 5 is entered to cbtain the concentration, I1, in per cent 
of helium by mass. 

The remaining step before pletiing of the concentration data is 
the determination of the centerline from: the rnaximam concentration 
point. This stepis carried out after the 4 p correction has been applied 
in order to determine the true maximum. In soe of the tests only 
enough points were taken vertically above and below the centerline to 
determine the true raximum for rate of decay versus axial downstream: 
distance. If a profile were somewhat sxewed, the centerline was 
deterrnined by equal spacing of half-maximum points. (This point will 
be covered in more detail in the later discussion. ) 

The helium: flow for sonic of the early tests was not exactly the 


0.0031 1lb/miin of the later tests. It was changed after the profile rung 
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to iruprove the exactness of the setting. <A scale factor of 1.15 was 
applied to the data of the early runs to refer all runs to the same helium 
flow. The linearity of such a correction for small ranges is shown by 
Figure 7, which is a plot of a few points taken at different helium flow 
rates with the probe fixed at an arbitrary location behind the model. 

This procedure would be in error if the concentrations measured were 
not in the linear portion of the millivolt versus per cent calibration curve, 
Figure 5. 

Finally, before proceeding to the results of the experiment, an 
estimate should be made of the errors to be expected in the concentration 
data produced by the equipment and procedures employed. The types 
and causes of the errors have been covered toa great extent in the above 
discussion. A more cornplete summary of them is included in Table IV, 
with an estimate of magnitudes wherever possible. This table may serve 
as a guide to improvement of the accuracy of the method. A fairly com-~ 
plete reference Concerning accuracy of thermal conductivity measurement 


is the book by paynest?. 


In 10. Summary of Test Parameter Variations and Test Objectives 


In view of the sin:plicity of the apparatus in the present exploratory 
investigation and the resultant long time intervals required in the 
measurement of individual wake samples, the number of parameters to 
be varied in the tests was neld to a minimum consistent with tae basic 
objectives previously cited. The most important parameter in the wake 
flow processes of interest is of course the Reynolds number. Variation 


of Reynolds number was achieved by holding the reservoir temperature 
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at To = 275°F and altering the tunnel reservoir pressure alone, 
according to Eq. lb. (Reynolds numbers of the tests are included in 
Table V.) The helium ejection rate was held at a constant low nominal 
value of . 0031 1b/min throughout the tests. This was ascertained to be 
the Zlow quantity desirable for proper tracing of the processes with 
minimum effect on the momentum of the basic flow. 

The investigation, with the above considerations, included the 
following tests: 

(1) Determination of the maximum concentrations behind the 
model for representative axial distances to the limits of probe travel 
(27 diameters) for supply pressures ey = 10, 35, 60, 85 psig, and all 
other conditions held constant. The maximum concentration would be 
determined by taking three or four samples near the tunnel centerline 
until the maximum is well defined. Tnese tests were required to 
determine how the maximum concentration decays with axial distance 
behind the model. 

(2) Measurement of flow samples for concentration at a sufficient 
number of points off the vertical centerline to determine a concentration 
profile at various downstream stations. These traverses for a number of 
such stations would then display conditions of spreading and regions of 
similarity in the flow. The combination of centerline maximum cone 
centration points taken and the profiles should then determine regions of 
turbulent flow. To restrict the amount of tunnel operating time, the profile 
traverses were limited to the extreme pressures considered, Po = 10 and 
35 psig. 


(3) Total pressure traverses were planned at a representative 
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downstream distance, x/D = 5, for the supply pressures P= 10, 35, 69, 
and &5 psig, with and without helium injection. Cornparison with data 
from the brass model would then show corrections that should be made 
for pressure differences, in utilization of that data for the porous model. 

(4) Total pressure traverses were to be taken off centerline to 
determine if the flow is two-dimensional at a represcatative station, 
x/Ds= 5 for SS 10 and 85 psig. 

(5) Monitoring of the flow at intervals during the concentration 
runs was also to be included, by the taking of sufficient total head 
traverses to see if conditions in the flow field remained consistent. 

These were considered the minimal tests that could be made to 
determine a satisfactory general description of the flow and diffusion 
processes affecting the injected helium. From sucha limited investi- 
gation, then, @ much more extensive survey could be planned, with 
improved equipment and technique for coverage in reduced tunnel 
operating time. 

Each test was designed to duplicate points of a previous test to 
check for repeatability of helium flow setting, model installation character- 
istics, conductivity apparatus variations, and possible alteration of the 
basic flow. Correlation of such repeated data would then establisha 
confidence level and accuracy check of the procedures. 

A typical data sheet for the recording of concentration test data 
is included as Table I. ‘The fixed parameters are included in the sheet, 
and the remaining blanks are explained by the discussion of the previous 


sections. 
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ll. RESULTS AND DISCUSSION 
Ill.1. General Discussion of the Plow 


At hypersonic speeds the flow field around a blunt body is domin- 
ated by the bow shock (Figures 8, 9). The static pressure on the front 
portion of a circular cylinder with ite axis perpendicular to gucha ctreai. 
is much greater than ambient pressure, and the pressure falls steadily 
along the surface. When the boundary layer is laminar, flow separation 
occurs sone distance aft of the position of maxinium thickness (Figures 
S, 9). The static pressure along the inclined free shear layer is very 
nearly constant and thie layer is nearly straight until Just upstrearn of 
the characteristic 'neck'' formed in the wake. Then the flow deflects 
and the pressure rises abruptly in a short distance, while the flow velocity 
in the wake cavity along the dividing ''zero" streamline is brought to rest. 
A second shock wave is produced by this flow deflection. This oblique 
wave is intercepted by the expansion fan venerated at the body surface 
and decays rapidly in the downstrean: direction. 

When the free shear layer ia laminar the width of the "neck" is of 
the same order as the boundary layer thickness at the separation point 
on the body. In the Reynolds number range of the present experiment, 
the cylinder boundary layer is laminar at separation, which is confirmed 
by the fact that the neck is extrem ely narrow in con parison to the body 
diameter. Rapid lateral diffusion of vorticity and heat takes place ag the 
fluid moves downstrean.. Of course, the gradiens of enthalpy and velocity 
in the ''external" flow which has traversed the bow shock are also bein# 


aroothed out by laminar diffusion, but the time scale of this process is 
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several orders of magnitude longer. 

li a foreign gas is ejected from the surface of the circular cylinder 
it will also be concentrated at the narrow ''neck'', provided the free shear 
layer is laminar. Downstrearn of the neck, the lateral diffusion of this 
tracer gas furnishes a measure of the efficiency of the laminar mixing 
process. At sufficiently high Keynolds numbers, transition to turbulent 
flow is observed to occur in this inner wake downstream of the neck. 

The subsequent lateral diffusion of the tracer gas in this case gives some 
insight into the rate of spreading of the turbulence into the ''external'' 

flow field. Of course the problem is complicated near the neck, but 
further downstream the lateral concentration profiles should exhibit a 
certain similarity. Consider a point far enough downstream of the neck so 
that uvcu<<u,. Then the Oseen approximation can be employed, and 


the diffusion equation to this approximation is as follows: 


pu, (8K/dx) = (8/2y p 2, ge) (9) 


where D,, = D>» the binary diffusion coefficient for laminar flow, and 
D,. = De » the turbulent eddy diffusion coefficient for turbulent flow. For 
laminar ilow me absorbs the density into the Howarth-Dorodnitzyn variable 


by letting 


Y -( (p/p) dy 


Oo 


By utilizing the approximation 
2 
pip D,>) = constant = Cp, (Dj), 


as analogous to the usual py = constant for the velocity defect solution, 


Eq. (9) is reduced to the form: 








C{(D, 5} om 
ane eee (10) 
e oY 


In addition there is the condition for the conservetion of the tracer gas 
Bas, 


namely, 
ca = \ pu xX dy = constant, (lla) 
90 
or co 
m = Pe a. ) IX GY = constant ° (11b) 


“Oo 
The solution to Eq. (10) with the condition given by Xas. (lia) and (11b) 


is the welleknown Gaussian distribution 


f. 








K = = exp (- apogee) 12) 
Ss wt a y> a nw Fe 7 
l2‘e @ = 
where 
A= A, (rh/p 2) , (A, = a constant) ; 
Te. , 
a y° 
K = A exp ( - ————————————— }) . (13) 
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The conclusions indicated by this first order laminar solution are that 
WK / Yz »givenp,,m 3 
(b) Since (D2). Yo Vales) ; 
en 1/ YP. , given x, m ; 


y* 


(c) K/K = exp { - Ae ela Ga ) ; 
max Ly - Ue 


indicating the Gaussian shape of the concentration profile. 


Now for turbulent flow the situation is more complicated because 
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the turbulent eddy diffusion coefficient is not independent of axial distance. 
Roughly speaking, De = Bau b-- » wheredau is the velocity difference 

across the turbulent inner wake, B = a constant, and b.- is the breadth 

of the inner wake. As the turbulent fluid penetrates into the "external flow", 
D,.. increases. Nevertheless a first approximation for the turbulent case 


is obtained by regarding D.. asa constant. Making another first order 


LT 
approximation, let p =p, a mean value in Eq. (9). Then the lateral 
concentration profile for the turbulent case is similar to the larninar 
solution, except that (Di. is replaced by Da. , ¥Y is replaced by y, 
and C= 1 in Eq. (13). Now 


he BIC) v 
where (C5). is the drag coefficient associated with the momentum defect 
of the turbulent inner wake and Dis the cylinder ciameter. Then for 


turbulent flow; 
exp ( gee) - (4) 


D ¥x/D 
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The conclusions drawn from the first order ''solution'' of the turbulent 


case are that 


(a) a as l/ Y= , given Pos» mM, ip) ; 


(b) K_47 WP,» givenx/D ,m, C, 


W 


me 


(c) The profile of K/K is Gaussian. 


max 
Before examining the experimental cata in the light of the 


preceding rough similarity considerations, certain flow conditions 





discovered during the tests must first be discussed. 
It. 2. Tunnel Flow Conditions and Variations 


During the course of the pressure traverse tests made with the 
brass model and the concentration tests with the porcus model, certain 
discrepancies began to appear. (Both models were mounted in the same 
tunnel position, 13.2 inches from the nozzle throat.) Discrepancics 
were first seen in the pressure data taken by McCarthy, prior to the 
concentration runs. (See Table VL) In particular, it was noted that 
reproducibility of the pressure traces was not good, and that square-shaped 
and non-symmetrical wakes were being obtained. Actually, Gaussian 
distributions of wake diffusion quantities are to be expected, according 
to the discussion in Section Iil.!. The discrepancies in the pressure 
traverses were attributed in pari to oil droplet formation on the brass 
model. {it diffused over the surface of the porous model.) However, 
it was not until a quantity of pressure and concentration data was taken 
that the discrepancy affected the data seriously enough to be certain that 
it, in fact, existed. The chronological history of the flow can be followed 
by referring to Tables VI and VII and the corresponding figures (Figures 
10-21), which are the complete presentation of the experimental data 
obtained. In particular, Figure 13 shows the effect of the flow change 
in spreading of helium mass; Figure 14 shows the effect on axial decay; 
and Figures 16-18 are pressure traverses which show the variations in 
the momentum defect of the wake core. Figure 19 demonstrates that the 
flow is identical for the two models. Figures 18 and 19 demonstrate that 


the effect of helium flow is that of a slight lowering of the wake total 





pressure, by 5 per cent or less at Py = 85, and widening of the shocks 
and wake core at ae 10 psig. Density effects of the light gas are 
undoubtedly the factors which caused these differences in total head 
raverses. 

When all of the acove figures are reviewed, one comes to the 
conclusion that tunnel flow perturbation is certainly evident. A few 
pressure traverses were taken without the model. These traverses 
and some calibration history of the tunnel revealed thai two small waves 
originating in the nozzle throat region crossed the tunnel at almost the 
exact position of the model. Their symmetry had made them impossible 
to detect in the early data. Their effect was to cause intermittent 
variation in the pressure of the wake and widening of all the wake total. 
pressure "'buckets''. This effect can be verified by noting that the 
pressure traverses taken with the probe off center actually correspond 
closely to those taken with the model rnoved rearward 2.6 in. out of the 
perturbation region. (See Figures léc, 18b, 18c, and 19b.) These con- 
clusions were confirmed by intermittently taking concentration data and 
pressure traverses with the model at the rear position until the data of 
concentration runs 10 - 15 was obtained without another shift of the flow. 
The runs taken with a stable flow then serve as a means of correlating 
the previous concentration profile data and cbtaining some useful 
information from it. This procedure is justified only because of the 
precautions taken in all concentration tests, namely those of repeating 
sample points from preceding runs to check reproducibility. The small 
variations in the duplicated points which were noted at the time of taking 
the profile data could now be understood. Originally, these were attri- 


buted to errors in helium flow setting or to instrument errors. 
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lt was considered impractical in this exploratory study to repeat 
the daia of the first runs for the present investigation, after moving the 
model to the stable position. Rather, it was decided to utilize the data 
in its present form to extract as much information as possible, and then 
to repeat and expand the experiment in another investigation, with im- 
proved equipment based on the knowledge gained to date. Actually the 
effect of this compromise on the remaining discussion is not as great 


as might at first be expected. 
lil. 3. Continuity Check of Helium Flow 


By integrating the diffusion equation to suitable limits of the 


concentration profile, one obtains 


+57. +6. +5, 

iS 9K 
| (2/ex(paniay + | @/eytpuwray = | (8/eyWeD,, Fray 
~ 84 “51, -5., 


Since the concentration and concentration gradient in the y- direction are 


both zero at the wake edges, the integral equation is reduced to 


57 
es 
(8/8x) | puk dy = Q ; 


x 
where the order of differentiation and integration has been interchanged. 


This result can only be valid if 


Ore 
j pukxk dy = constant = rays unit length ‘ (15) 
“5. 
which expresses the fact that the mass of helium injected from the cylinder 
into the stream must equal the integral of the puk profile taken at any 


station, x, and integrated over the span. Thus if one assumes two- 


dimensionality of the flow, the integral 
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b Si 
{ { puk dy dz = Tie » mass/sec (16) 
° “5. 


would equal the helium flow quantity metered into the cylinder. The 
integral provides a means of checking whether or not the concentration 
measurement is accurate, providing the product pu can be determined. 
The method of computation for pu feliae from the equation of state and 


appropriate modifications: 


pu = p(Ma/RT) = pM Y(¥m/7ROT) Ss (7) 
where Ym and Ri, are the ratio of specific heats and gas constant of the 
mixture of helium and air. If the assumption is made that the flow is 
iso-energetic (with possible error not greater than approximately 10 per 


cent), the local temperature can be replaced by the well-known relation, 


T> 
kn 
1+ ———— oi 


The Mach number may be determined from the pressure data taken by 
McCarthy, which corresponds to the concentration profiles taken before 
large shifts in the tunnel flow. The method of obtaining the Mach number 
from this data is given in Section I. 8, for Y= Y¥_. Although the static 
pressure is given only for the wake center, McCarthy has found by taking 
representative traverses that static pressure is very nearly constant over 
wake cross~gsections between inner shocks downstreain of the wake ''neck", 
To provide greater accuracy in the puk calculation, one should 
apply a correction to the static pressure for the helium defect in the wake. 
(See Figures 18 and 19.) A measurement of the static pressure with and 


without helium flow should ordinarily have been made. It is reasonaoble, 
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however, that the velocity is not changed anc that the static pressure, p, 


me 
5 PRS a 


because of the density effect of the helium. The effect in thig case was 


is changed by the same ratio as the total pressure, P 


so slight that no correction to p was made. The Mach number does not 
change in any case, with the assumption that the velocity change is very 
small, 
An accurate computation of the puk quantity also requires that 
the parameters YmandR_ be replaced by their true values, rather than 
m 
. . , y = Z 4 Ole 
by approximating them with Ja=1.4and R, = 1,716 ft. “/sec” ~R for air, 
because the parameters oi the two binary constituents are quite different in 


the case of helium and air. (R = 12, 438 ft. “eee °R and Vise = 1. 667) 


Hi 


The cquation for the gas constant is 


Rn = y KR, = KR, tU-K) RR (19a) 


beg ad 


where K = concentration of helium inair, by mass. The equation tor the 





ratio of specific heais is cy 
, : me 
(l-K)c_ +Ke (1-K)F, + KE ie om 
Y = Fa PHe . (20a) 
Yi “K) CL AC. ae ° 
: ie l=K + K —seee 
3 ve 


When the numerical values for the gas parameters are substituted, the 


two equations become 


Rin 


1,716 + 107. 22% (19b) 


Yn, = 140 f 5. 251K (200) 


The value of puk obtained for each vertical distance, y, by the 
above procedure, is plotted asin Figure 23. The intecration by planimeter 


then gives the mass per unit span. The double integration is complete when 
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the mass/span is multiplied by the span length, b= 5 in. The two 
example curves of Figure 22 came from the profile data of Figures 10 
and 11 and the corresponding pressure and Mach number profiles of 
Figures 20 and 2l. The non-symmetry of the puX profiles is to a certain 
extent the result of non=-symmetry in the total pressure profiles, but 
may also be caused by the non-uniform injection of helium through the 
cylinder. Four such integrations were Gb wEibLe with the available data 


and were completed with the result: 


P, = 85 x/D = 15 Th, = + 00279 Lb/min 
P, = 85 x/D = 20 Thy, = + 00300 lb/min 
P, = 10 x/D = 9 Misr, = - 00307 lb/min 
P, = 10 x/D = 15 rhs, = - 00281 lb/min 


These values are to be compared with the flowrator setting Mize = . 00305 
= .00010 lb/min. (The second value above was taken by doubling a half- 
concentration profile, after the flow had shifted slightly.) Even though 
the several approximations above were made, the values have less than 
1Q per cent error and 10 per cent spread. It could be argued that the 
results are somewhat fortuitous, but it is unlikely that gross errors 
could have cancelled each other for the four different cases at two tunnel 
reservoir pressures. One is led to the fairly safe conclusion that: 

(a) The iso-energetic temperature assumed is within about 10 
per cent of the true value. 

(b) The static and total pressure measurements are dependable. 

(c) The concentration measurements are also to be trusted 


within the accuracy predicted by Table IV. 


(dq) The tunnel instability was not serious in this case because 





the pressure and concentration profiles corresponded. 

(ce) Even better accuracy should be obtained if the instrumentation 
is improved and the model is located in the rear position. 

({) Two dimensional flow is, for all practical purposes, a 
certainty. (This conclusion is supported by the relatively good agreement 
between the forward off-center pressure traverses and those taken wita 


the model in the rear position. ) 


li. 4. Discussion of the Diffusion and Similarity Properties 


The values of all maxinium concentrations taken at the wake 
centerline are plotted in Figures 14 and15. The plot of Figure i4 was 
included to show the effects of the flow shift oa axial decay at P, = 85 psig, 
when the model was in the forward position. The points plotted in Figure 
15, however, are those taken with stable flow conditions when the model 
was moved to the rear position. Unfortunately this shift in model position 
reduced the probe travel behind the model centerline to 27 diameters 
instead of the previous 36, but by this expedient a complete and cone 
sistent set of axial decay points at four Be values was obtained. The 
flow was checked before, during, and after the runs by total pressure 
traverses with the probe at 5 diameters. The stable traverses of Figures 
16 - 19 were consistently duplicated. Having established a fairly high 
degree of confidence in the data of Figure 15 by the continuity check of 
the previous section and the evidence of stable flow, this plot rnay be 
discussed in the light of previous experiment and theory. 

When the word similarity is used in the following discussion it 
should be construed to mean similarity or approximate similarity in the 


general sense; that is, a scaling of the flow according to certain laws 
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relating the coordinates and diffusion parameters, such that all cross- 
section profiles taken at downstream positions have the sarne structure 
and shape and follow the same decay rule. Townsend attached this 
meaning to the term " selfepreserving''* , and reserved the terzn 
"similarity" for the more specific Reynolds number similarity in 
turbulent flow, requiring that the ''processes which determine the main 
structure of turbulent motion are substantially independent of viscosity". 
This essentially applies to all turbulent free boundary flows in the con- 
tinuum regimes at high velocity. It does not mean, however, that the 
boundary conditions of the flow cannot be influenced by viscosity, as 
indeed they must in the case of the cylinder wake, where the drag 


coefficient is a function of Reynolds number. 


UI. 4.1. Downstream Decay of Maximum Concentration 


The theory shows for the present case that the maximum helium 
concentration should be distributed down the tunnel centerline in inverse 
proportion to the square root of x, the axial coordinate. Further 
ar a5 should be true of the flow in either the laminar or turbulent 
case. The plot of Figure 15 indicates that this condition is true in only 
the Po = 10 psig case for the region beyond x/D = 4. A constant slope 
for the other reservoir pressures is indicated, however, for regions 
sufficiently far downstream. in the physical case, it is realized that the 
model axis is not necessarily the origin of the similar flow, and that in - 
order to compare the slope of the decay curve with theory, one must 


determine the virtual origin of the similar flow by a suitable auxiliary 


plot. This procedure is accomplished for the present data by plotting 
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the values 1/K against the coordinate, x, ona rectangular plot. 


max 
(See example, Figure 23.) Though scatter is enlarged by this procedure, 
a straight line through the points is evidently a good approximation to 
their distribution, showing that the 57 condition is satisfactorily met, 
with virtual origin as indicated by the x-intercept. The resulting 
correction to the plot of Figure 15 for all four reservoir pressures is 
included in Figure 24. [It can be seen that the rnaximurn concentration 
does at least closcly approximate the decay condition given by the theory. 
It should be noted that a few points in Figure 15 fall below the 
normal scatter of the data. A close check of the experimental data sheets 
revealed that the points n.arked in the figure were ‘alcen at the end of a 
day of operation with nigh tunnel make-upair. It is quite likely that the 
tunnel drying equipment could not remove sufficient quantities of water 
vapor in these cases, and the conductivity of the tunnel air was lowered 
as compared to the reference sample taken at the beginning of the run. 
Nith the relative conductivity of water vapor as given in Table I, about 
1] per cent by mass would have been sufficient to cause the observed 
deviations. This is thena likely reason for the appearance of the profile 
taken at 24 diameters in the P, = 10 psig case of the early data (Figure 11) 
which was also taken at the end of a long test period with high make-up 
air. The use of high make-up air for other than short runs therefore 
appears questionable. Another expedient would be to collect a reference 


sample for each tunnel sample, but this would increase the time of 


taking a reading and the possibility of error in the procedure. 


-_ 
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Iii. 4.2. Variation of Maximum Concentration with 


Tunnel Pressure 


The separation of all the curves in Figure 24, for the various 
reservoir pressures, does not appear to be a function of P. alone, but 
is dependent on some other condition of the flow. Since the theory calls 
for K~ 1/ YP. proportionality at a given x/D in the laminar case, one 
immediately checks to see if the condition is satisfied in the present data. 
This is done by checking the slope of the auxilliary plot made to determine 
the virtual origin, xX» OF by checking the ratio between straight-line 
portions of the curves in Figure 24. The condition is satisfied for the 


case of Po = 10 versus Po = 35 psig, with 


where 14.4 psi is the average barometric pressure of the laboratory 


elevation. Itis not satisfied for 


YP 17 Po } — rps 17 Pt : 
where the values for the absolute square root ratio are 1. 23 and 1, 16, 
respectively. Furthermore, the flow does not satisfy the predicted 
cendition for the fully turbulent case in either of these comparisons, 


in that the iG 2 values are not scaled directly by the pressure ratios, 
P ,(60)/P.(35) and P 9(85)/P,(60) 


A possible conclusion indicated by these facts is that the flow at 
Fe = 60 psig is ina region of transition from laminar to turbulent flow 
and that the flow at P. = 85 psig is fully turbulent. The flows at P, = 10 


and 35 psig are of sufficiently low Reynolds number that they are laminar. 
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These statements arc supported by the relatec prsssurc traverses which 


were taken intermittenily with the concentration suns to check flow 


S 


stability (Figures 16 - 18, runs 10 - 13). Thc wake portion of the pressure 


traverse in the a = 60 psig case is narrower aac deeper than that of the 
Po = S85 trace. This trend is opposite to that between the Ps = 10 and 


Eo 


35 psig cases, where the Po = 10 trace becomes shallower at the 
wake centex on the nondimensicnal scale, and remains at about the same 
width. Schlieren photograpns and hot-wire studies made by A. Demetriades 
also support these conclusions. Transition ai 3 - 4 diameters in the 
Po = 65 case is indicated, and at 4 - 5 diameters i. tne P, = 690 psig case. 
in both cases it appears ihat the mechanism of transition cormin.ences in 
the 2 diameter region,corresponding to the position of the necked-down 
portion of the wake. (See Figure 15.) The other possible cxplanation for 
the lack of predicted scale between the i a turbulent = curves, is that 
the drag aifects the scaling in a maaner which cannot be determined by 
the present data. (The appearance of CH in the boundary conditioa of 

Re W 
the turbulent case was snuown in Section fi.1.) The present catia is 
inadequate in that momentum thickness coulc not be properly cormputea 


because of the effect of the tunnel pertuz bation. 


lJ.4.3. Transverse Spread of Concentration 


The similarity of coacentration profile shapes may be checked by 


plotting the normalized concentrations K/K versus the normalized y 


max 
coordinate a /y » where Pe ig the mean of the positive and negative 
absolute values and Y,, i8 the same mean taken at the half-maximum 


point of concentration. The curves are thus fitted at two points and the 


scatter pattern of the remaiing points checked for similarity. The plots 
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of Figures 25 and 26 are the result of sucha procedure. It can be seen 
that, except for the x/D = 3 case at P= 10, the curves indicate a 
similarity as to shape of the spreading concentration profile. The 
exception noted in the x/D = 3 curve is to be expected, since it is too 
close to the origin of the similar region and an anomaly appears in the 
data of Figure 11]. The profile points taken at = 35 after a slight 
shift of the flow, plot equally well with the profiles at 9 and 15 diameters 
before the flow shift. Having thus established at least reasonable shape 


L 
similarity and the decay with x 4, one may plot K/K, oe versus 





and obtain one curve for all profiles. This is done for the P, = &5 case 
in Figure 27, with the appropriate virtual origin for the flow shift. In 
addition, a mean, 7 , and standard deviation, © , are computed from 
the curve faired through the points, and the two parameters are used to 
compute a Gaussian distribution for comparison with the faired curve. 
The mean is just the width of the half-profile base and the standard 
deviation is given by 


coe 7° (6/Kna sd 5 
" Ss wena (i 


for discreet values of the abscissa. The chi-squared test of statistical 
theory could be applied to determine the percentage confidence in the 
normality of the distribution. It was not done for this case because the 
tunnel perturbation effect would alter the meaning of the result. 


Reasonable certainty of agreement with the theory, however, is indicated. 
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In spite of the tunnel flow diificulties, it can be stated in summary 
that the flow is approximately similar for x/D = 9 throughout the Reynolds 
number range investigated. The related incompressible literature 
indicates that much greater downstream distance is required for the 
flow to achieve complete pinnae Rity®. One possible reason for this 
difference is that the proper transverse length scale for the hypersonic 
"inner wake'' is actually much smaller than the body diameter, so that 
a distance of 9 - 10 body diameters corresponds toa distance of many 


inner wake momentum thicknesses. 
Ill. 4.4. Non-Similar Regions of the Flow 


Very little quantitative evaluation of non-similar flow regions 
has been presented in the literature. This is understandable in view of 
the fact that these are regions of rapid adjustment in which nonlinear, 
fluctuating effects take place in the transition from one state of natural 
flow to another. In the wake of the hypersonic cylinder, such regions 
are to be found in the narrow portion of the wake and in the near vacuum 
at the rear of the model. A typical profile observed in the wedge shaped 
Cavity between the model and the neck of the wake is that of Figure 26. 
The trapezoidal shape is evident until the concentration profile achieves 
a triangular distribution at the neck, after which the statistical effects 
of diffusion alter the shape toward that of the error curve. The base 
and point of the triangle are last to achieve the characteristic normal 
distribution in the similar flow downstream. This effect is also observed 
in the incompressible coaxial jet profiles of Forstall and Shapiro! ’. The 


trapezoidal region corresponds to their jet potential core. 
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aT}. Sugg estions for Improvement of the Experimental Method and 


Scope of the Wake Investigation 


ly. 8. 1. Improvement of the Apparatus 


Several suggestions have previously been made concerning the 
apparatus, but a few additional points could be made, after reviewing 


the data obtained. 
Il.5.1.a. Model Construction 


The data reveals that porosity of the model was adequate, but 
improved uniformity would possibly remove some of the unsymmetric 
behavior of the profiles. A few variations in the flow were noted which 
are not accounted for by the tunnel perturbation. It is possible that 
these resulted from model surface condition effects on the drag. The 
surface conditions are affected by model tool marks and cther defects 
and the oil problem in the tunnel. These factors can be reduced toa 
ininimum if the oil problem can be solved and the same model and model 
orientation is used for each test. The concentration measurements are 


very sensitive to drag effects on the wake pressure. 
Ili.5.1.b. Probe Inetallation 


The probe is restricted in the present installation to about 27 
diameters of travel, after moving the model rearward. The traversing 
mechanism should be modified to give about 50 - 60 diameters of travel 
downstream. It should be provided with a positive locking system to 
prevent shift during measurements, anda micrometer scale. It should 


be small enough in cross-sectional area to provide easy starting of the 
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flow. (This was a problem in the present investigation, with a probe 


mechanism of about 5/8 inch circular cross-section. ) 


lI.5.1.c. Helium Flow Metering 


The metering of helium. flow in the present investigation was 
satisfactory only after it was realized that niodei hack pressure caused 
variances in the rm.eter reading for the same mass flow, if the sr.all 
meter was used. The large meter provided the necessary stability to 
differences in model porosity, but it was necessary to read the scale 
very accurately. The reading could be improved by the installation of 
a light and magnifying glass on the meter, with the eee of a means 
of sighting the level of the ball through a reticle or gunsight arrangement 
to insure consistency of eye level. With this arrangement the helium flow 
could be turned off and on at will and readjusted quickly, which would 
reduce the usage of helium during periods of apparatus adjustment and 


check-out. 


Il. 5.1.d. Sampling Apparatus and Procecure 


There are two major difficulties with the present sampling system 
that should be corrected before extensive testing is resumed: 

(1) The times required to clear the previous sample and bring 
in a new one are unnecessarily long. 

(2) The accuracy of the pressure correctionis not satisfactory 
for very low sample densities. 

For the first problem, there are two possible solutions: The 
first has been previously tnentioned as the pumping of each sample to 


barometric pressure. This method has attendant difficulties with regard 





56 


to possible contamination of the sample. [In addition it requires a valve 
sequencing and pumping method which may become complex. It does, 
however, also eliminate the second problem. The second method is to 
follow the procedure utilized in the present investigation, measuring the 
pressure defect after a single Compression, but to use larger compression 
ratio and provide automatic sequencing wherever possible. The automatic 
sequencing is desirable for precision and alleviation of the natural 
tendency to allow extra time for each step of the procedure. It could be 
accomplished by a suitable system of relays and solenoid-operated 
valves. The time would be even more reduced if a sensitive recorder 

is used in place of the potentiometer and the pen dropped for each sample 
point, in the relay sequencing. The proper time constants could be set 

in the relays with a few preliminary measurements at each reservoir 
pressure. The system and valves would have to be leak checked, of 
course, at frequent intervals. A rough diagram ofa suitable systern is 
included as an Appendix. 

The second problem, of measuring the pressure accurately, can 
be handled by reading total pressure intermittently with sample readings 
on the same recorder plot, if the sample is not compressed. Compression 
is desirable, however, at the very low tunnel densities, to prevent the 
pressure correction from becoming as large as the sample reading. If 
a transducer were installed on the sample volume the pressure could 
still be read intermittently on the recorder, and the total pressure at 
the probe could also be plotted as a third trace. 

The proper combination of the above ideas should provide equipment 


which would operate to take sample points automatically after each 
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movement of the probe. The Appemlix is an example of one such possible 
combination. Sample points could be read in times of less than one 
minute if equipment of this type could be developed, instead of the 
several minutes average per point required in the tests of the present 
experiment. 

In addition to these problems, the water vapor contamination is 
of some additional concern. Either high tunnel make up air should be 
utilized and a reference sample gathered from the tunnel circuit, or 
low make up air with the single reference sample as taken in the present 


tests. 


UWI.5.l.e. Instrument Specifications 


The potentiometer of the present experiment was of Leeds and 
Northrup manufacture and was sensitive to four decimal places in 
millivolt units. It was well suited to the tests. Three decimal place 
accuracy is sufficient, but stability of zero is important. A recorder 
utilized in place of the potentiometer would need to meet the same 
specifications. 

The ammeter utilized was of the standard laboratory variety, 
accurate to about one-tenth milli-amp. An accuracy of 0.01 milli-amp 
is desirable for maintaining constant current through the cell bridge. 

The conductivity cell itself was fairly accurate and stable, but 
it should be cleaned at intervals and checked frequently for possible 
drift of zero or potential bias of the setting made to read the initial 
calibration curve. Known samples are indispensible to the proper 
employment of the conductivity cell, and gooc mixing equipment and 


proceedure should be employed for their preparation. The current and 
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voltage used in the cell bridge should be as high as possible for sensitivity 


but low enough to prevent fusing of the paired resistance wires. 


Lise. Suggestions for Additional Coverage of the 
Cylinder Wake 


Any new equipment designed for concentration studies in the ake 
should be checked against the results of Figure 15, which is the plot of 
points taken with the model in stable flow. Following check-out of the 
equipment, the profile data of this investigation should be repeated in 
the stable flow. The new data should include K nax points taken for 
various tunnel pressures in the turbulent Reynolds number regions. 
Extensive total and static pressure coverage of the same Po and x/D 
points should be made to determine exchange coefficients. Extensive 
coverage of the necked region of the wake should be obtained to determine 
how initial conditions might affect the later sirnilar flow. Comparison 
of the normalized profiles with error distributions and modified 
formulas from the empirical theory should be made. Variation of the 
helium flow is desirable to see what effects may result in the diffusion 
and pressure profiles. Another diffusing gas should be employed to 
determine the effects in this regard. A heavy gas for instance might be 
useful in determining the relative diffusion of mass and momentum as 
compared to the light gas of the present experiment. An extensive program 
of wake coverage as outlined will provide a good background for studies of 


all types in wakes and boundary layers, where binary gas components 


may be detected by the thermal conductivity method. 








DY 
TV. CONGLUSIONS 


The several important conclusions which may be derived from 
the present investigation are as follows: 

(1) The thermal conductivity method is an important and useful 
one for diffusion studies in the hypersonic flow field, where one binary 
constituent is to be detected in another. The method is adequate in 
determining mass diffusion properties of a flow witha light gas injected 
as tracer or as coolant for the model. 

(2) The problem of low densities in the samples taken from 
hypersonic flows can be controlled by suitable calibration. The procedure 
can be improved over that of the cxperiment with a more sensitive 
measurement of the sample pressures. 

(3) Other improvement of the model, probe, and instrumentation 
will increase the accuracy of the method employed in this case to within 
one or two per cent error for concentrations on the order of one-half 
per cent. 

(4) The method is adequate to determine turbulent processes of 
the flow and to detect transition, if the data is properly reduced and 
plotted to bring out the similarity properties. 

(5) The wake of the hypersonic cylinder is formed from the 
boundary layer gas which is retained initially in shear layers at very 
low pressure behind the model, in trapezoidal shaped profiles of mags. 
The gases are compressed into a very narrow wake at a distance of 
about two diameters downsiream of the model, after wiich the lateral 
diffusion rapidly approaches the similarity regime. 


(6) The similar pattern is set up in either laminar or turbulent 
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flow at between 5 - 9 diameters behind the model. This result is in 
contrast to the much greater distances required in low speed flow. 
(7) In the laminar case, the theoretically predicted scaling of the 


maximum concentration inversely as the square root of the supply 


pressure was satisfied in the ratio of tne P, 10 and P, = 35 psig 


readings. In the higher pressure cases (Po = 60, 85), K ~ aid not 


ma 
scale with ee as an approximate consideration indicates. Either the 
dependence on the pressure is more involved in the turbulent case, or 
the fully turbulent flow was noi yet reached at the highest Reynolds 
number of the experiment. 

(8) Diffusion profiles were affected somewhat by a tunnel flow 
perturbation, but are still observed to be typically Gaussian, as theory 


predicts, a few diameters downstream of the narrowed portion of the 


wake. 
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SUMMARY OF GAS PROPERTIES 


Molecular Ratio of Temperature Per Cent 
Weight Conductivity Coefficient of in Dry 
te Air at Conduction Air at 
0°C 0 - 100°C, Sea Level 
per C By Méess 


Air 1.00 - 00265 


Carbon Dioxide . 605 . 00527 
5.97 . 00256 
7.145 . 00265 
-996 . 00264 
1.013 , 00303 
685 .00311 
.725 (46°C) . 00540 





Column References: 


(2) Handbook of Chemistry and Physics, Ist edition 


(3), (4) Daynes, Gas Analysis by Measurement of Thermal Conductivity, 
Table 1, Cambridge University Press, 1933 


(5) Aeronautical Handbook, Pratt and Whitney Aircraft Company, 
1957 
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TABLE H 


CONCENTRATION TEST DATA SHEET 


Run No. Date 


Type of Run 


(1) Tunnel Conditions ad ie ini ! obe ample ern! 
y, in. |} Ap, Cm Hg! mv 


Me 5. & eg = 275°R 


Fs = psig 


(2) Model installation Check 


Vac. pump vs. charnber 


Vac. pump vs. model 





Difference, (cm Hg) 





(3) Room Conditions 


Bar. Press. Temp. Time 
cm Hg Cc 








ess, P= He. 
(4) Instrument Check 
zero bridge tare riniv 


e en ieee iene 
Sample correction ratio, tor 
bridge bias 


(5) Helium Flow Data 


Hie on at : psig 

‘He off at : psig 

Flowrator setting 

Metering pressure cm Hg 
j | O 

Metering temperature Cc 








Tunnel Conc., mv Time 
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TABLE IV 


ERROR ESTIMATE 





1. 

helium low sensitivity 
flow of flowmeter 
setting 

2. 

Zero unstable T. C. 
shift cell, ammeter, 
and 

potentiometer 
a 


pressure| low sensitivity 
calibras | of pressure 
tion scale 

error 


3) ts«(4 


Amount Range of 
Error 


23.5% all readings 
of emf 
(mv) 


= 05 all readings 


mv 


(5) 
Possible 
Remedy 


more sensitive 
meter 


stable sensitive 
instruments; 


rnonitor current 
closely; keep 
cell clean 


use microe 
manometer on 
pump scale; 
measure Ap to 
.1 mm Hg 


4, Miscellaneous small effects having less than 1G aggregate error 


for average concentration readings of the tests: 


probe setting 


pee at 


oil in flow 


temperature difference e. vibration 


{. heliurn impurity (from lines) 


model differences g. tunnel conditions off by 


- .2épsiande F 


5. Intermittent errors (large ones can usually be detected): 


a. water vapor in flow caused by excess make-up air with 
insufficient drying 
b. model installation or heliurm line leakage after metering 


c. large contimination in conductivity cell 


dad. bad leaks in sampling system 


Total possible error for a sample of 1% helium concentration at pressure 
ratios ap/Ph . <.9 is approximately 67%. 
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TABLE V 


TEST REYNOLDS NUMBERS 


F = constant 


Re/in x Vole *e 





* Based on atmospheric pressure of 14. 4 psi for the laboratory 
elevation, and cylinder diameter, D= 9.30 inch 
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TABLE Vil 


SUMMARY OF CONCENTRATION DATA 


Protile 
11/30/59 Profiles 
12/3/59 Profile r . 00279 lb/min 


12/22/59 Profile 
12/23/59 Profiles r . 00307 lb/min 
00231 Lb/min 


1/29/60 K 


‘max Flow shiftea slightly 


since runs 1-5 


2/2/60 


Ka pts 
2/4/60 K 
2 


pte 
max 
Profile 2 


max 





ie a 00300 1lb/min 
fo! x/iD= 20 


2/29/60 pts 


1 
**mase 


w 
4, 
Profile 15 
4 
3 





Runs 10-]5; 
Basic flow shifted 


10 14,15 3/18/60 85 XK 


Me 3/21/69 10 Rag DE 2. ? again since runs 6-9. 
12 15 3/22/60 10 K nay Pts 5,20,27,1 Model position now 
35 Sh. cu 2.6" aft of position 
13.15 3/23/60 35 KO epee 8, 27) 972 EL EN 
60 2,5,9, 20 


i Mp6 values refer to continuity check. See Section II. 4. 


This table is continued on the following page. 
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TABLE Vile- CONTINUED 


(5) (7) 


P.. Type of D Remarks 
psig un 





14 15 3/24/60 60 K pte 5, 3,7, 


ae 15, 27 
Profile ] 
14,15 85 Kets Grp09. o,8h2, 
sa 1.5, 18, 27 


15 14,15 3/25/60 35 K ax Pte 3.1 
85 pee 
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/— 8 feed- through connectors wired 


in full Wheatstone bridge. Fiber- 


glass insulating sleeve,!2 leads. 


For chromatography to 300°C. 
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Gow-Mac Instrument Co. 
Madison, N.J. 

Mode! 9285 with wiring 
mod. 9193(Te- I!) 
Tungsten filaments 


Dimensions of cell: 
Bix 2x 2 5/6 
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DIAGRAM OF WIRING FOR THERMAL CONDUCTIVITY 
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FIG.§5 MASS CONCENTRATION OF He IN No VS. BRIOGE 
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FIG.6 CALIBRATION OF CONDUCTIVITY CELL FOR PRESSURE 


DIFFERENCE BETWEEN SAMPLE AND REFERENCE 
CHAMBERS 








ie 





0.00410 
0.00305 4 
0.00220; — x —_—++-_--- ee | ey oe | a 
| | 
on | 
Si Ordinates correspond to | | 
tr min flowmeter scale readings | | 
| | 
0.40 0250 0.60 070 0.80 O.90 


Concentration, K,m.v. at arbitrary position in wake 
of cylinder injecting helium 


FIG. 7 LINEARITY OF CONCENTRATION WITH A SMALL 
CHANGE IN FLOW OF HELIUM 








HiGeeo =- SCHLIEREN ee OF POROUS CIRCULAR 
CYLINDER, M2 3-8 42t Pere ee  L = 275°F 
(extraneous pattern caused by oil in flow; model 
shape affected by chips in glass) 





FIG- 9 -- SCHLIEREN PHOTOGRAPH OF POnRODOUS CIBRGU Lavi 
CYLINDER, M=5.8 at P. = ag psi¥ TO = Opera 





Jl 


6 | 
Concentration , K, “% 
14 
2 Runs 1-3 > M=5.8 


Por 85 psig, To= 275°F 


m = 0.003! Ib/min 
He 









AD 
Bee: 
A 5 
oOo 9g 
Vv iS 


FIG. 10 HELIUM CONCENTRATION PROFILES IN THE WAKE “OF 
A POROUS “CRFCULAR CY UIN@ER 





Bz 





Concentration, 










K, %o 
Re 
| | 
| P ir Runs 4, 9 
oo nr eee eee LOD, Suing 
1 ome ocr 
M = 5.8 
| Me? 0.003! Ib/mig, 
! ‘ x/p : 
| /\ a O35 
oe es ND a 
| { A 9 
Galo 
fp -V 
| S VY 24 (water vapor lowered 
| A a concentration ) 
\/ 
| Y Approximate effect 
7 \ \ / of water vapor 
. : oe at x/D=24 
P] O — eee eae 
—| O | 
YD 
oer | HELIUM CONCENTRATION PROFILES IN THE WAKE 


OF A POROUS CIRCULAR CYLINDER 





So 


Concentration , 7 


Kanes | 












? 


: 
Runs 8,9 | 


Rog Bors 
To= Af Deak 
(Slight tunnel flow shift since fig.!0) 
M=5.8 
m1. 70.0031 1 / min 


«/ 
O 20 (half profile) 
A 19 


re 





¥/ 1.0 2.0 
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APPENDIX 


THERMAL CONDUCTIVITY APPARATUS AND PROCEDURE 
FOR HYPERSONIC WIND TUNNEL EMPLOYMENT 


Refer to Figure A-1 


I. Preliminary, with probe at desired distance behind model: 

1, Check all lines, volumes and valves tor leakage. 

ec. Check concentration measurement on known sample. 

3. Start tunnel and stabilize. (Can be done along with 1 and 2.) 

4, Make zero concentration check with air in sarm:ple and 
reference cavities. Zero the concentration scale on the 
recorder. Tura helium on. 

5. Prepare to set pressure sensitivity by manually opening 
valve a, c, andd with m, b, andeclosed. Measure total 
pressure on the micromanometer manually. Closed. Zero 
the pressure scale on the recorder. 

6. Set total pressure sensitivity in recorder ordinate scale; 
find tunnel centerline by probe movement vertically and set 
probe position potentiometer to read correctly on abscissa. 

W. Measurement of Total Pressure and Concentration: 

A. Automatic Operation, where sample is at pitot pressure 
greater than 1 cm Hg absolute: 

1. Set probe vertical distance, valves a, c, e, and f 


open, and band d closed. Turn on ''washing machine'"’ 
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type relay system, which follows the following 
sequence; 

{2} Closes vaive m and opens n to clear cell - 
relay holds for sufficient time, then opens 
m to allow next sample to wash through - 
relay again holds and then closes m. This 
operation leaves the cell clean and at low 
vacuum for drawing of the tunnel sample. 

(b) Relay operates after sufficient time to achieve 
cell vacuum and closes n, then immediately 
opens im. (With l cm Hg or more differential, 
a tunnel sample is drawn into the cell in 
satisfactory time for automatic operation. ) 
After sample is drawn, relays switch to 
transducer circuit ana drop recorder pen 
momentarily for total pressure, then switching 
to cell bridge circuit, drop pen for concentration. 

(c) Relay cycle is complete. Move probe and repeat 
for other points. This procedure obtains both 
total pressure and concentration. 

B. Manual Operation, where sample is at pitot pressure less 
than 1 cm Hg absolute: 
The pressure correction in this range is too great for 
direct concentration measurement of the sample at pitot 
pressure. The sample must be compressed. The hand- 


operated pump is adequate in that very few readings are 
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taken at these low pressureg with the preseat tunnel 
densities. The rmercury pump itself may be designed asa 
micromancmeter, in which case it may be used for the 
reading of cell pressures and will be more convenient for 

the manual operation than the transducer. (For a description 
of the hand operation, see Sections IJ. 7.1. and III. 5. of the 
basic report. ) 

The measurement of total pressure along with compressed 
sarnple pressure requires additional steps. Usually, however, 
it will not be necessary to obtain total pressure at these low 
sample pressures, as they occur only in non-similar regions 
near the model. 

If possible, the cell bridge current should be turned off 
during each interval when the cell is completely evacuated. 


(See Section I. 7.1. of the basic report.) 
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